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We consider here a full study of stellar dynamics from the brane-world point of view in the case of 
constant density and of a polytropic fluid. We start our study cataloguing the minimal requirements 
to obtain a compact object with a Schwarszchild exterior, highlighting the low and high energy limit, 
the boundary conditions, and the appropriate behavior of Weyl contributions inside and outside of 
the star. Under the previous requirements we show an extensive study of stellar behavior, starting 
with stars of constant density and its extended cases with the presence of nonlocal contributions. 

Finally, we focus our attention to more realistic stars with a polytropic equation of state, specially 
in the case of white dwarfs, and study their static configurations numerically. One of the main 
results is that the inclusion of the Weyl functions from braneworld models allow the existence of 
more compact configurations than within General Relativity. 

PACS numbers: 04.50.-h,04.40.Dg 


I. INTRODUCTION 

Stellar astrophysics is one of the most characteristic 
topics studied by General Relativity (GR), which has 
helped to describe the dynamic and evolution of stars 
with unprecedented success {TJ. In addition, the matter 
inside a star may be in some cases in extreme conditions 
generating complicated high energy phenomena, princi¬ 
pally in white dwarfs, neutron stars, and others, and then 
a complete description of the stellar properties requires 
the introduction of a particular equation of state (EoS) 
like in the case of polytropes[2], or even Bose-Einstein 
Condensates (BEC)|S]. 

Another interesting possibility in recent times is to 
consider alternative theories of gravity and to look for 
their particular signatures in stellar models, specially for 
some of the extreme situations mentioned above. For 
instance, the authors in[4] considered the corrections in¬ 
duced by a Galileon Lagrangian in stars of constant den¬ 
sity. Another example is given by the so-called models of 
Brane-Worlds (see [3 [ja for a good review) whose main 
characteristic is the existence of branes (four dimensional 
manifolds) embedded in a five dimensional bulk|7b This 
particular geometry allows a natural extension of Ein¬ 
stein’s equations 0: and introduces new degrees of free¬ 
dom through quadratic terms of the energy momentum 
tensor, the non-local Weyl terms, and other fields that 
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could live in the bulk. This framework has been used for 
stars with a constant density in{9], and also for polytropic 
matter with a given relationship between the quantities 
arising from the non-local Weyl terms in[T(]:. It has also 
been shown that the exterior solutions of these brane- 
stars is not the Schwarzschild onefS, IT] , and then the 
Weyl fluids in the exterior of the stars can have a non- 
negligible influence in the internal pressure and compact¬ 
ness of stellar objects. More recently, the conditions for 
stellar stability in brane-stars were revisited in[I2| for a 
set of hypotheses called the minimal setup , which are 
consistent with a Schwarzschild exterior. Also see^3] f° r 
a study on the gravitational collapse of brane stars. 

With the previous background, this paper is dedicated 
to the study of the stellar equations of motion that arise 
from the formalism of Brane-World theory, and the role 
of the Weyl functions in the regular behavior of a stellar 
distribution. It is important to remark that our main 
objective is to consider models of stars as realistic as 
possible, and for this reason we will follow conventional 
wisdom in this regard: a Schwarzschild exterior, and reg¬ 
ularity of all functions involved. Based on these premises, 
we perform numerical studies of the so-called extended 
GM solution with constant density, and of a polytropic 
fluid. 

The organization of the paper is as follows. In 
Sec. [II] we describe the equations of stellar dynamics 
with branes, emphasizing the high and low energy limits, 
boundary conditions, and the role played by the Weyl 
functions in providing consistent and regular solutions. 
Subsequently in Sec. |III[ we study the case of constant 
density and the extended GM solution. Also, in Sec. [TV| 
we study polytropic brane-stars. Finally, in Sec. [V] we 
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give some conclusions and remarks. 

II. STELLAR DYNAMICS WITH BRANES 


where a prime indicates derivative with respect to r. We 
have also defined V = 6 U /, AC = 4V /, and A(r ) = 
[1 —2G N A4(r)/r] _1 , whereas p e ff and p e ff are explicitly 
given by: 


Let us start by writing the equations of motion for 
an embedded brane in a five dimensional bulk using the 
Randall-Sundrum II model[7]- We first assume that the 
Einstein equations are the gravitational equations of mo¬ 
tion of the 5-dimensional Universe, 


Gab + A^gAB — K (5 )Tab ■ 


(1) 


Following an appropriate computation, the modified 
4dim Einstein’s equation can be written as|B}|F 

G T C fiu T T ^F^ , (2) 


where k ( 4 ) and K( 5 ) are respectively the four and five- 
dimensional coupling constants, which are related one to 
each other in the form: k 2 4 ^ = 8i tGn = ) A/6, A is 

defined as the brane tension, and Gjy is Newton’s con¬ 
stant. For purposes of simplicity, we will not consider 
bulk matter, which translates into F = 0, and dis¬ 
card the presence of the four-dimensional cosmological 
constant, A( 4 ) = 0, as we do not expect it to have any 
important effect at astrophysical scales (for a recent dis¬ 
cussion about it see[H]). Additionally, we will neglect 
any nonlocal energy flux, which is allowed by the static 
spherically symmetric solutions we will study below [5]. 

In the case of a perfect fluid, the energy-momentum 
T^ v tensor, the quadratic energy-momentum tensor II M1/ , 
and the Weyl can be written as: 


— pu^u v -f- ph^ , 
n Mi , = ^p[pu^u v + (p + 2p)h ^), 


S/w — 


.2 ’ 
'(4)' 


Uu u ii u + Vr u r u + 


U-V 


(3a) 

(3b) 

, (3c) 


where p = p(r) and p = p{r) are respectively, the pres¬ 
sure and density of the stellar matter of interest, U is the 
nonlocal energy density, V is the nonlocal anisotropic 
stress scalar, u a is the fluid four-velocity, that also satis¬ 
fies the condition g^u^iF = —1, and h^ = g^ v + 
is orthogonal to u ll . Under the assumption of spherical 
symmetry, the metric can be written as: 


ds 2 = — B{r)dt 2 + A(r)dr 2 + r 2 (d9 2 + sin 2 9dip 2 ). (4) 
The equations of motion for stellar structure then are[(T, 

H2J: 

M' = 4:Trr 2 p e ff , (5a) 

P 1 = ~\^{P + p), (5b) 

V' + 3AC' = -^-(2V + 3A0 - 9 A/" - 3 (p + p)p ’, (5c) 
Jd r 

B' _ 2.G n f 4nr 3 p eff + M \ 

B r 2 \l-2G N M/r J ’ { ’ 


p,// =p(i + f) + fx + 

/’•ii = A 1 + jx) + T 


V 

3A 


AC 
A ’ 


(6a) 

(6b) 


Now, we are in position of analyze the following impor¬ 
tant points. 


A. Numerical analysis 


In order to have a numerical solution of the equations 
of motion, we choose the following dimensionless vari¬ 
ables: 


x = y 7 G n M/R ( r/R ), p = p/(peff), P = p/{p e ff){7 a ) 

A = A/ (p e ff), V = v/(p eff ) 2 , AC = A r/(p ef f ) 2 , (7b) 


for which Eqs. © read 
M' = x 2 p eff 


P =-2 

X* 


V' + 3A V = —„ 


3 ( X 3 p e ff + M 
1 — 6 Ad / x 
X 3 p e ff+M 
1 — &M/x 


(p + p ) > 

(2V + 3AC) 


(8a) 

(8b) 


—AC — 3(p + p)p ', 

X 


(8c) 


where now a prime denotes derivatives with respect to 
x, the mean effective density is just given by (p e ff) = 
3M/4irR 3 , and also 


Peff = P 
Peff = P 



AC 

A ' 


(9a) 

(9b) 


Note that the ratio p/\ is invariant under the change 
of variables, p/X = p/ A, and then we will omit the bar 
whenever the ratio appears in the equations of motion. 


B. High and low energy limits 


There are two very clear limiting expressions of the 
equations of motion in terms of normalized brane ratio 
A, as the latter represents the energy ratio of the brane 
tension with respect to the mean energy density of the 
compact star of interest, see Eq. (7b). It is usually as¬ 


sumed that the brane corrections are measured in terms 
of the absolute value of the brane tension A, but in our 
study we find the brane ratio A = A /(p e ff) to have a more 
meaningful character for compact objects in general. 
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Under this line of reasoning, we first present the low 
energy limit of the equations of motion, represented by 
the operation A —>• oo, under which Eqs. ^ become the 
usual Tolman-Oppenheimer-Volkoff (TOV) equations of 
GRUJE]: 


M! = x 2 p , 



x 3 p + M \ 
1 — 6 M. jx) 


{P + P ), 


(10a) 

(10b) 


where the effective pressure and density are directly rep¬ 
resented by their normalized physical values. 

We have called it the low energy limit because we 
are assuming that the mean density of the star is much 
lower than the brane tension A, so that any brane cor¬ 
rections in the equations of motion are highly suppressed 
by the brane energy scale. We cannot say here whether 
the brane tension is at a very energy scale, or it is just 
that the star density is not high enough. In strict sense, 
Eq. ( |8c| ) can still be considered for the integration of the 
Weyl functions, but their values will not make any dif¬ 
ference in the final integration of the physical variables. 

There is also the high energy limit of the equations 
of motion represented by A —> 0, for which the effective 
density and pressure read 

p 2 V 

PefS ~ 2A + I ’ (lla) 

P.// = ^(2f + « + b +j- (Hb) 


We can see that there is an overall factor of A in the 


above expressions (11), which will also appear as such in 
Eqs. ([8]) in the high energy limit. The brane ratio can 
then be absorbed in the equations of motion by means 
of the following change of variables: M. —> A4A 1 / 2 , and 
x —¥ x\ 1/2 , and then we finally find the equations of 
motion for the high energy limit: 


M'=^(p 2 + 2V) , 


r x 3 (pp + p 2 1 2 + V/3 + N) + M 


P =-2 


x* 

(p+p)-. 


1 — 6 M/x 


(12a) 

x 

(12b) 

-tv , okp 6 f x 3 (pp + p 2 /2 + V/3+Af) +M\ w 

V +A/V 1 -6M/x J 

(2V + 3J\f) — —J\T — 3(p + p)p'. (12c) 

x 


C. Boundary conditions 


The change of variables 0 is very appropriate to ex¬ 
plore the solutions of the TOV equations ([8|, as all phys¬ 
ical quantities involved are normalized in terms of two 
important observables in stellar astrophysics, which are 
the mass M and the radius R of the star. Furthermore, 
these two parameters appear in the single combination 
GjvM/i? that represents the compactness of the star. For 
instance, the interior range of the new radial variable is 
x — [0, ^GnM/R], which means that the surface of the 
star is located at x(R) = X = -\JGnM / R. Also, the new 
mass function changes to: 


M(x) 


1 /G,vM\ 3/2 M(r) 

3 [ R J M 


(13) 


and then the total mass is M(X) = (1/3)(GjvM /i?) 3 / 2 . 
In other words, the compactness of the star will deter¬ 
mine the mass and size of the numerical solutions. 

The equations of motion will be integrated from the 
center up to the surface of the star defined by the con¬ 
dition p(X) = 0; the latter only refers to the physical 
pressure, and we will take it as a reasonable physical as¬ 
sumption even though it is not necessarily required in the 
case of brane stars. Finally, at the center of the star we 
will also assume that A4 —► 0 as x —» 0, so that there is 
not a discontinuity of the different quantities in the cen¬ 
ter of the star, and the central value of the pressure (or 
any other related quantity) will be set as a free parameter 
that will characterize the numerical solutions. 

Even though we will not consider exterior solutions, we 
must anyway take into account the information provided 
by the Israel-Darmois (ID) matching condition, which for 
the case under study can be written asjU: 

(3/2 )p 2 (X) + V~(X) + 3AT~(X) = V+pO + 3AT+(A), 

(14) 

where the superscript — (+) denotes the interior (exte¬ 
rior) values of the different quantities at the surface of 
the star, and we also assumed that p(x > X) = 0. 

A desirable property we want in our solutions is a 
Schwarzschild exterior , which can be easily accomplished 
under the boundary conditions V+(X) = 0 = A/" + (A'), as 
for them the simplest solution that arises from Eq. |8c| 
is the trivial one: V(x > X) = 0 = J\f(x > A'). Thus, for 
the purposes of this paper, we will refer hereafter to the 
restricted ID matching condition given by: 


In contrast to the TOV equations of GR in (10), we 


shall call this case the high energy limit because the mean 
density is much larger than the brane tension, even if 
we cannot say whether this is because the brane tension 
attains a very small energy value, or it is just that the 
star has such a large density that the latter surpasses the 
energy scale of the brane tension. 


(3/2)p 2 (A) + V-(A)+3AT-(X) =0. (15) 

For completeness, we just note that the exterior solu¬ 
tions of the metric functions are given by the well known 
expressions B{r) = A^ 1 (r) = 1 — 2 GnM/x. In addition, 
it can be shown from Eq. (5d) that the interior solution 


of the lapse function, in terms of the normalized vari- 
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ables ([7]), is given by: 


B (x) 


1 - 2G n M/R 


= exp 


r X 


dx 


6 


X 3 p ef f 


M 


x 2 \ 1 — 6 M/x 


( 16 ) 

and then we will not solve it explicitly in any of the cases 
presented below. 

The numerical recipe described above will be applied to 
different cases and configurations in Secs. |rnland|TVl The 
results that will be obtained will have a universal char¬ 
acter, as they will not depend upon the particular values 
of the mass and radius of a given star, but they will rep¬ 
resent general classes of stars according to their common 
compactness GnM/R. This will allow us to reach wide 
general conclusions about the physical properties of the 
different configurations by means of numerical methods. 


D. Weyl functions 

It must be noticed that the interior solutions cannot 
evade the presence of the Weyl terms even if the exterior 
solution is Schwarzschilcl. For example, let us put by 
hand that Af(x) = 0. If the density is constant p(X) ^ 0, 
the ID matching condition (15) implies that V~(A) = 


— (3/2 )p 2 (X), and then the full solution must be[9!: 
V(x<X) = ~(3/2)p\l+p/p) 4 


(17) 


The full consequences of this nonlocal energy density V 
are explored in Sec. |IIIB| below. 

Even the condition of a Schwarzschild exterior together 
with p(X) = 0 do not directly imply that the Weyl func¬ 
tions must vanish in the interior, as it can be shown P"2 : 
that in such a case Eq. (8c I must have the following so¬ 
lution: 


V(x < X) = 


B 2 (x) 


r x 

/ B\p 

J X 


+ p)p' dx , (18) 


which accomplishes the boundary condition V~(A') = 0 
and is also regular at x = 0. This is particularly impor¬ 
tant for all cases in which the density is not constant, as 
we shall see below for the polytropes in Sec. m 

In the opposite case when V(x) = 0 and the density 
is constant, the ID matching condition ( fl5| ) implies that 
at the surface of the star: Af~ {X) = — (T/2)p 2 (A), and 


then Eq. (8c I integrates into 


A f{x <X) = ~- 


1 f X 


{p + pf 


(19) 


There is though a non-divergent interior solution of Af 
if we drop the condition of constant density, which is: 

N(x < X) = * J Bx 3 (p+p)p'dx. (20) 


But the ID matching condition (15) now indicates that at 


the surface of the star we must have p(A) p 0, or either 
give up the Schwarzschild exterior. In consequence, the 
only interior solution of the nonlocal anisotropic stress 
under the conditions of a Schwarzschild exterior, and 
non-constant density with p(X) = 0, which are the con¬ 
ditions we expect to have in realistic stars, is the trivial 
one: Af(x) = 0 (see also [T2] ). 

There are other possibilities that have been explored in 
the specialized literature, like for instance a relationship 
between the Weyl functions in the form J\f = <rV, where 
0 is a constant parameter [10j. Clearly, the solutions 
and (18) are special cases for which o = 0. In the general 


case, Eq. (8c) can be written as: 


(1 +3(7)1/ = -^-(2 + 3cr)V- -aV -3{p + p)p', (21) 
B x 

as long as <7 p —1/3. If the density is constant, then 
there is a solution which is similar to Eq. (19): 


V = C (p + p) 


2(2+30-) x ~9a 


1 1/(1+3<t) 


( 22 ) 


where C is an integration constant that could be deter¬ 
mined with the help of the ID matching condition (151. 


However, this solution is not appropriate for the interior 
of the star because, as it happened too for Eq. (19), it 
diverges in the center of the star (x = 0). 

We can also consider the case of non-constant density, 
in which case we find a similar solution to Eq. (20): 

-| 1/(1+3<t) 


Af(x < A') = B(x) 


\-2(2+3a) x ~9a 


jg 2(2+3cr) a .9tr 


-1 1/(1+3ct ) 


(p + p)p'ddf23) 


Even though this solution is well behaved in the interior 
of the star, it needs a non-trivial boundary condition at 
the surface as dictated by the ID matching condition (151, 


and for that we require either to have p(X) p 0, or to 
give up the Schwarzschild exterior. 

We see that the imposition of a Schwarzschild exterior 
has strong consequences for the interior solutions of the 
Weyl functions, mostly because it is difficult in general to 
find for them a well behaved interior solution. The most 
problematic case is that of the anisotropic stress function 
Af, and for this reason we will not take it into account as 
part of the brane gravitational corrections, but assume 
that the latter are only given by the quadratic corrections 
of the density p 2 and the nonlocal energy density V. 


Needless to say, this solution diverges at the center of 
the star and cannot be considered as a useful interior 
solution. That is, in the case of constant density there is 
not a regular interior solution with the only presence of 
the nonlocal anisotropic stress Af. 


III. THE CASE OF CONSTANT DENSITY 

One of the simplest possibilities of star modes is that 
of constant density p, which can be solved under different 
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gravitational schemes. In this section we will work out 
such a case within the braneworld scheme and explain the 
additional physical and boundary conditions that may be 
needed in order to reach well posed numerical solutions. 


A. The case of the Germani-Maartens solution of 
brane stars 


To start with we consider here the GM interior solu¬ 
tion, which was thoroughly studied in[S], and that does 
not take into account corrections induced by the Weyl 
terms: V = 0 = A 7. The modified TOV equations 
are given again by Eqs. (5a) and ( |5b| ) with the follow¬ 
ing identifications: p e // = p( 1 + p/ 2A), and J? e // = 
p( 1 + p/X) + P“/2 A. 


Because the density is constant, we find, in terms of 
the variables in ([7]), that (p e ff) = p{ 1 + p/ 2A) and then 
p e ff = 1. Likewise, we find that p = (1 + p/2A) _1 , and 
then its value is directly determined by the ratio p/X. 
Notice that p < 1, and that we recover p = 1 in the GR 
limit p/X —>■ 0. The boundary conditions depend upon 
the compactness of the star G^M/R, as in the case of 
GR, but also upon the ratio p/X, as expected in brane 
models. The exact solution of the pressure function is[S]: 


V _ y/1 ~ 6 M/X - y/1 - %Mx 2 /X* ^ 

P y/1 - 6 Mx 2 /X 3 - 3C" V 1 -6M/X’ 


the effective density and pressure: 


( 1+ 2X, 

\ 3 _p 

>-2 P X 

(1 +p/p)\ 

(25a) 

(-0 

1 

IQ-Ioi 

+ 

P 2 P x^m*. 

(25b) 


Peff = P 


Here we have taken into account that the nonlocal en¬ 


ergy density is given by Eq. (17). As it can be seen in 


Eqs. (25), there are negative contributions in both the 


effective density and pressure originated from the pres¬ 
ence of the Weyl nonlocal energy, and the solutions now 
depend upon three separate parameters: the constant 
density p, the brane ratio p/X, and the compactness of 
the star GnM/R. 

As we expect to have p{x) > 0, then the effective den¬ 
sity must be an increasing function, p e fj(x) < p e //P0. 
which may even attain negative values at the interior 
points where the pressure is largest. Moreover, we also 
infer from this information that (p e ff) < Peff(X) = 
p( 1 — p/ 2A), and then we expect that in general p > 
(1 — p/2A)W In contrast to the GM case above, p cannot 
be given a fixed value beforehand and becomes a variable 
that must be adjusted appropriately so that the numer¬ 
ical solutions accomplish all boundary conditions. This 
time, however, the GR limit p = 1 is a lower bound as 
p/ A —* 0, which is an early indication that the known GR 
upper bound on the star compactness could in principle 
be surpassed by the new solutions. 

The equations of motion are more easily solved if we 
take the following change of variables: 


where ( = (1 + 2p/A)/(l + p/X). The brane ratio p/X 
lowers the maximum value of the compactness of the 
star, and the numerical solutions satisfies the analytic 
bound found from the exact GM solution: G^M/R = 
(1/2)(1 — ( 2 /9). The GR limit is obtained when p/X —> 
0: GnM/R < 4/9, whereas in the opposite direction 
p/X ->oowe obtain: GnM/R < 5/18. 

According to the discussion in Sec. |IID[ the GM solu¬ 
tion cannot be matched to a Schwarzschild exterior, and 
for that reason it is usually assumed that other exterior 
solutions with the presence of the Weyl function must be 
the correct ones for brane stars. However, it has been 
recently shown [12]. under very general conditions, that 
the GM solution plays also the role of being the limiting 
case of realistic stars when brane corrections are consid¬ 
ered, and then gives an upper bound in the compactness 
of stars with both brane corrections and a Schwarzschild 
exterior. 


B. The extended GM solution 


We will now review the interior brane solution with 
constant density, a Schwarzschild exterior, and a non¬ 
null Weyl term V. This cases was also briefly considered 
in[5], but lacks an analytical solution. The equations of 
motion are again (|8j) with the following expressions for 


xp~ 1/2 


M-^Mp x l 2 ,w = p/p, (26) 


where w is the EoS, and then the Eqs. (|8|) become: 


M'=x 2 


i + A_?qi 

2A 2 A v 


■ w 


V eff 


t 2 
Peff 


X 3 W ef f 


M 


x* \ 1 — QM/x 


(1 + w ), 


(27a) 

(27b) 




The only free parameter that appears explicitly in 
Eqs. ( f27| ) is the brane ratio p/X. Moreover, the outer 
boundary conditions, see for instance Eq. ( |13[) , must 
be adjusted to the values X = {G^Ad/R) 1 ' p 1 ^ 2 and 
M(X) = {GnM/RY^p 1 ! 2 /?,. 

Examples of the numerical solutions allowed by 
Eqs. ( f27| ) are shown in Fig. [I] for the brane ratios p/X = 
10 1 ,10 6 , where it is confirmed that there are numeri¬ 
cal solutions well beyond the GR limit of GnM/R < 4/9. 
We only considered cases in which the star has an over¬ 
all positive mass, for which it must also have a positive 
density at its surface. The latter can be translated into 
the condition p e ff(X) > 0, and then from Eq. (25a) we 
find the constraint p/X < 1. 

There are two main reasons for the surpass of the GR 
limit. The first one is that the extra free parameter p is 
only bounded from below, and then it is at our disposal to 
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find numerical solutions that can surpass the GR limit for 
any given value of the brane ratio p/X. Correspondingly, 
the second reason is that the effective density p e ff is 
an increasing function that can become as negative as 
necessary in the interior of the star. Actually, as far as 
the numerical experiments are concerned, the only true 
limit that could be found for the numerical solutions is 
the Schwarzschild one GatM/J? <1/2. 


n as 7 = (n + l)/n. For example, white dwarfs can be 
modeled by the polytropic index n = 3, and neutron stars 
by polytropes with an index in the range n = 0.5-ing. 

The equations of motion | 8 ]) can be simplified if we 
follow the usual recipe for polytropes and make the fol¬ 
lowing change of variable for the density: p = 0 ", where 
n is the polytropi c ind ex defined above. For the reasons 
explained in Sec. IID we set Af = 0. Eqs. ([ 8 ]) are then 
written in the form: 



0.15 


0.1 

0.05 

0 

pA=10~ 6 




-0.05 

- \ 

- 

2 -0.1 

\ 

- 

-0.15 

- 

- 

-0.2 

G n M/R 0.30 
-G n M/R =0.35 

- 

-0.25 

G n M/R 0.40 . 

G n M/R 0.45 


-0.3 

G n M/R =0.48 

G n M/R =0.49 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 


x 


r / 8 n \ 

V 

0 ( 1 - 1 —= 

+ 

L V 2a; 

A 


3 I'x 3 p e ff + M\ (l+KO) 
x 2 \ 1 — 6 M/x J K(n + 1 ) 


12 / X 3 p e ff + AA 
2 \ 1 — 6 M/x ) 


(28a) 

(28b) 


9 / x 3 p eff + M \ 

x 2 \ 1 — 6 M/x J 


n (l + K0y q2u-i 
K(n + 1) 


,(28c) 


and the effective pressure (|9b| now reads 


Peff = K0 n+1 



d 2n y 
2A + 3A ' 


(29) 


It must be stressed out that in our case the density pa¬ 
rameter 9 gives an indication of the values of the density 
p with respect to the mean value of the effective density 
( peff ), given by the dimensionless density p, in contrast 
to the standard case in which the value of reference is the 
density at the center of the star p( 0 ). 

It can also be shown that, again like in the standard 
case of polytropes, the polytropic coefficient K is a re¬ 
dundant constant and can be hidden in the equations of 
motion. If we further consider the following change of 
variables: 


x -I xK n/2 , 6 8K - 1 , M -t MK n / 2 , (30a) 

Peff -> PeffK~ n , A -> A K~ n , V -> VK~ 2n ,(30b) 


then Eqs. (28) simply read 


FIG. 1. The profile of the integrated mass M(x) correspond¬ 
ing to the extended GM solution with ratios p/X= 10 _1 (Top) 
and p/X = 10~ 6 (Bottom). We can see that the extended GM 
solution allows the existence of stars with a compactness be¬ 
yond the GR limit but below the extreme Schwarzshild limit 
GnM/R < 1/2. Notice that one reason for that is that the 
mass function can acquire negative values in the interior of 
the star for the most compact cases. See the text for more 
details. 


IV. POLYTROPIC BRANE STARS 


In this section we study brane stars with a polytropic 
fluid and an EoS in the form p(r) = Kp^(r). Here, K is 
the polytropic constant, and 7 is the polytropic exponent, 
which can be written in terms of the polytropic index 


r,„ ( 8 n \ 

V 

8 n 1 + ^ 

+ T 

L V 2 Aj 

A 


3^ / x 3 p e ff + Ad A (1 + 9) 
x 2 \ 1 — 6 M/x ) (n + 1) ’ 


12 / X 3 p e ff + Ad A 
2 V 1 — GM/x ) 


9 / x 3 p e ff + M \ n(l + 6>) 2 2n _ 1 
x 2 \ 1 — 6 M/x J (n + 1 ) 


(31a) 

(31b) 


(31c) 


where 


Peff=0 n+1 


8 2 


V 


1 + ^ I+ ^A + 3A' 


(32) 


Our main interest are the numerical solutions of stars 
with a finite size, as determined by the boundary con¬ 
dition p(X) = 0, which in the case of the polytropes 
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translates into p(X) = 0, and from this into 0(X) = 0. 
In order to avoid any singularities in the equations of mo¬ 
tion at the surface of the star, in particular for Eq. (31c I, 
we must constraint the values of the polytropic index 
in the range n > 1/2. Needless to say, such a con¬ 
straint does not exist either in the case of non-relativistic 
(Newtonian) or relativistic (GR) polytropes. It must 
be noticed as well that the boundary conditions should 
also be adjusted so that X = (G^M / R) 1 / 2 K~ n / 2 and 
M(X) = (1/3 ){G N M/Rf/ 2 K- n / 2 . 

We now include brane corrections with the contribu¬ 
tion of one of the Weyl terms, with the boundary condi¬ 
tion V(_Y) = 0, so that the ID matching condition (14) 
allows a Schwarzschild exterior for the polytrope and dic¬ 
tates that the interior solution for the nonlocal energy 


density is given by Eq. (18). As discussed in Sec. IIB 


the brane terms must contribute to the effective density 
and pressure inside the star, which means that we can¬ 
not have in this case a counterpart of the GM solution, 
unless the Schwarzschild condition were waived. 

As in the cases studied in Sec. Ill we will integrate 
inwards the equations of motion under the same bound¬ 
ary conditions presented in Sec. m with the central 
value of 0(0) being a free parameter that will help us to 
classify the numerical solutions. The most compact star 
will be given by the maximum in the plot of the com¬ 
pactness as a function of the central density: GnM/R 
vs 0(0). The value of the compactness will be read off 
from the outermost points of the numerical solution as 
G^M/R = 3M(X)/X, whereas the polytropic coeffi¬ 
cient can be calculated from: K = [SA^Y)] 1 / 3 / X. 

This time we have to give explicit values to the brane 
tension A and the polytropic index n. For the latter, we 
consider in the following sections the case of white dwarfs 
with n = 3, whereas the brane tension will remain free 
to label the different brane star solutions. 


A. Numerical solutions 



FIG. 2. Numerical solutions of the int erio r profile of the 
mass M of polytropic brane stars, see Eqs. (311, with A = 10 2 . 
We can observe for GnM/R > 0.180 that the mass becomes 
negative, like in the constant density case shown in Fig. [l] 
This is due to the contribution of the Weyl function V, see 
text for more details. 


the effective density at the center becomes negative for 
GnM/R > 0.175. It is clear that the geometric term V 
contributes notoriously for large values of the compact¬ 
ness in the high energy limit. 

Thus, at least for the range of compactness we nu¬ 
merically explored, the contribution of the Weyl tensor 
affects the internal configurations of the stars in such 
a way that there is no maximum for the compactness 
that can be reached, except for the Schwarzschild bound 
GnM/R < 0.5. This can be seen in Fig. [4j where we 
note that for A = 10 6 ,10 5 the curves reach a maximum 
value just as in the case of polytropic stars in GR. How¬ 
ever, as A decreases it is possible to find stellar configu¬ 
rations with a larger compactness beyond the standard 
GR bound. 


Solutions for the high energy limit with A = 10 2 al¬ 
lowed by Eqs. (311 are shown in Figs. [2] and [3j In par¬ 
ticular, the interior mass profile M(x) is shown in Fig. [2] 
for a range of compactness: GnM/R = 0.170 — 0.230, in 
which the main feature we can observe is a change in sign 
close to the center of the star for G^M/R > 0.180. This 
behavior is due to the contribution of the Weyl function 
V in Eq. (31a). Note that the same behavior occurs in 
the constant density case in Section |III B[ so this type 
of effect from the Weyl function is present also in more 
realistic stars. 

On the other hand, numerical solutions with A = 10 2 in 
Fig- 0 show that for low compactness the effective den¬ 
sity has the expected decreasing behavior as we move 
outwards from the center of the star. However, as the 
compactness increases the maximum value of the den¬ 
sity is displaced from the center, and then the density 
profile is not just a decreasing function. Also note that 


V. CONCLUSIONS AND REMARKS 

In this paper we studied the equilibrium configura¬ 
tions of stars with gravitational corrections in braneworld 
models, and provided numerical solutions when neces¬ 
sary. For that we considered the high and low energy 
limits of the equations of motion to show the threshold 
between GR and braneworlds, and explored the appro¬ 
priate boundary conditions to obtain general conclusions 
about the physical properties of the different stellar con¬ 
figurations. 

Our analysis took into account the corresponding Weyl 
functions which provide non-local terms in the pressure 
and density, and which can have noticeable effects in di¬ 
verse features of a star. This study allows us to relinquish 
the nonlocal anisotropic stress under the conditions of a 
Schwarzschild exterior and non-constant density, which 
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FIG. 3. (Top) The effective energy density p e // as a function 
of the radial coordinate x, which results from the numerical 
solution of Eqs. (311 for low compactness and A = 10 2 . It 
can be seen that, as the compactness increases, the maximum 
of the effective density is displaced from the center of the 
star. (Bottom) The interior profiles of the effective density 
p e ff for high compactness. The effect of the Weyl term V is 
sufficiently large to change the sign of the effective density at 
the center; actually, p e ff(x = 0) < 0 for G^M/R > 0.175. 


are conditions rightly expected for a real star. 

As initial test, we revisited the case of constant density, 
corresponding to the GM solution, but later studied the 
so-called extended GM solution, for which exist stars with 
a compactness beyond the standard GR bound. This is 
due mainly to the existence of the non-local terms which 
provoke the appearance of negative values of the effective 
density and mass in the interior of the star. 

Finally, we considered the case of a white dwarf star 
modeled with a polytropic EoS and index n = 3. In sim¬ 
ilarity with the extended GM case, our results proved 
the existence of dwarf stars with a compactness beyond 
the GR limit, because of the presence of non-local terms. 
Also, the compactness of dwarf stars is not bounded as 
the brane tension tends to zero, which correspond to the 
high energy limit, while in the low energy limit we recov¬ 


FIG. 4. The compactness GnM/R as a function of the 
central value 0(0) for the polytropic configurations obtained 
from Eqs. ( |31[ ). It can be seen that the compactness of the 
polytropic star with n = 3 is not bounded as A — > 0. For 
A = 10 6 ,10 5 the curves coincides with that of polytropes in 
GR. 


ered the classical compactness reported in the literature 
for the case of GR. All these results are in agreement 
with the study in[T2], where it was shown that one of the 
main assumptions for the existence of an upper bound in 
the compactness of a star was that the effective energy 
density in the interior should be a decreasing function. 

The results presented were based on a clear method¬ 
ology that make the equations of motion of braneworlds 
more tractable in numerical terms that in other analy¬ 
ses in the literature. As noted, the presence of brane 
corrections modifies in a notorious way the compactness, 
mass, and other physical characteristics in stellar dynam¬ 
ics, even under the assumption of a Schwarzschild exte¬ 
rior. Extended studies along the lines suggested in this 
paper can be used to constrain the value of the brane 
tension using observational data provided by stellar dy¬ 
namics, and with that to find evidence for the presence 
of extra dimensions. As a final note, we cannot say if all 
the configurations found would be gravitationally stable, 
but it is very likely that those with a negative effective 
density in the interior may not be able to prevent the col¬ 
lapse into configurations well within the general bound 
found in|12j. This is work in progress that will be re¬ 
ported elsewhere. 
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